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ABSTRACT 


In the measurement of microwave circuits and devices, the data reduction of two- 
port measurement data obtained at the external terminals of microstrip-coupled devices, 
with correction for measurenient error, can be obtained by de-embedment procedures 


with use of the HP 8409B vector network analyzer and a Fortran program. 
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I. INTRODUCTION 
A. BACKGROUND 


Microwave circuit elements, such as semiconductor devices, cannot be connected 
directly to the terminals of a microwave network analyzer for measurement of scattering 
Or impedance coefficients but must be coupled through a chain of hardware elements 
comprising connector-adapters, microstrip sections, and bondwire connections. The 
basic problem under consideration is the removal of the effects of intervening hardware 
from measurements made at the external terminals of an embedded device or circuit 
component to obtain the intrinsic device characteristics. Accurate measurement data 
on a device is essential for successful circuit design. 

Earlier work on device de-embedment has been coupled with undesirable restrictions 
such as the assumption of physical symmetry of hardware, inadequate provision for re- 
alistic measurement standards, or neglect of statistical treatment of data for improve- 
ment of accuracv. Consequently, a need exists for a practical treatment of this 
measurement problem which can be implemented in typical laboratory circumstances, 
with acceptable accuracy of results. 

In order to measure the scattering parameters of a device, the device is often 
mounted in a fixture. This fixture has a set of scattering parameters associated with it 
so that the parameters measured at the terminals of the fixture are the scattering pa- 
rameters of the device transformed by the embedding network of the fixture. 

Thus, to obtain the scattering parameters of the device itself, the device must be 
“de-embedded” from the fixture. Techniques for using a network analyzer (N.A.) to de- 
embed the device from the test fixture have been discussed earlier [Refs. 1,2,10,11]. 

As long as the device interface connection is coaxial line or waveguide, de- 
embedding is relatively straightforward because standards are available. But if the device 
is constructed without a coaxial or waveguide interface, as in the case of microstrip or 
coplanar line systems, de-embedding becomes more difficult. 

Extracting the device parameters from the overall measured parameters is strongly 
dependent on the ability to characterize the fixturing. Figure | is the block diagram and 
signal flow graph of a typical N.A.- fixturing setup. The problem ts to find the best 
characterization networks A and B so as to remove their contributions, yielding the de- 


vice parameters, S. 
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Figure 1. Block Diagram and Signal Flow Graph of a Typical Fixturing Set-up 
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B. SCOPE 

The scope of the proposed research will be to develop a working procedure for car- 
rving out the removal of embedded device scattering matrices from network analyzer 
data. The procedure developed will be supported by experimental measurements on 


tvpical microwave circuits fabricated on planar substrates. 


Il. CHARACTERIZATION OF THE FIXTURES 


A. CHARACTERIZATION CHOICES 

Several papers present characterization techniques that are in use at various labo- 
ratories. There are four main choices of N.A.-fixturing characterization as shown Fig. 
OReie al 


CHARACTERIZATION CHOICES 


Veto Method 2 


(1). Short ckt 
fitted to a 


(2). 4 Open ckt 
model! 


Method 3 Method 4 


(1). Through - More than 3 
(2). Short known 


(3). Delay ref'coeff’ of 


(Speciale) each port 





Figure 2. The Four Main Choices of N.A.-Fixturing Characterization. 


Method (1) combines a variety of measurements, including the transmission between 
fixtures connected back-to-back and reflection measurements on the fixtures terminated 
in microstrip stubs. In this way, data is assembled to enable the scattering coefficients 
of the embedded device to be converted to the device scattering matrix. Method (1) 
therefore contains features of both methods (3) and (4) below, and represents also an 
early example of the use of microstrip stubs as known reference standards. 

In method (2), the measured S-parameters are fit to an assumed circuit model to 


calculate the device parameters. 


Method (3) has the advantage that only three simple test pieces in the device me- 
dium are required. More details are given in Chapter V, part B. 

From the above choices, method (4) was selected for investigation. In this method, 
three or more known reflection coefficients are measured at each port. The reference 
load reflection coefficients are derived from independently measured or calculated non- 
zero values. Reference 2 used a chip capacitor, wire bonded to microstrip on alumina, 
with capacitance measured at low frequency (LF). The overall reflection coefficient was 
then measured versus frequency. To obtain multiple calibration measurements, a piece 
was scribed off the chip, and the measurement cycle repeated twice more at each port. 

In other cases, the reference load was a high-Q varactor chip, again measured as a 
function of three values of reverse bias at LF. The chip was then wire-bonded to the 
microstrip conductor and the overall reflection coefficients of the fixture-chip combina- 
tion versus frequency were measured at each value of reverse bias. 

In the present work, the reference load is a section of the microstrip line itself. Three 
different lengths of microstrip line yield the three different reflection coefficients for each 


port by separate measurements. 


B. THE TEST FIXTURE 

The S-parameters of the test fixture can be determined by measuring a set of known 
impedances in place of the embedded device. This process is known as “unterminating.” 
Once the fixture’s S-parameters have been determined, the device is measured in the 
calibrated test fixture, and the device parameters are “de-embedded” from the data to 
give an accurate picture of device performance. 

The signal flow diagram and block diagram in Fig. 3 define the input and output 
characteristics of the test fixture. 

Three sets of S-parameters must be determined when an unknown device is charac- 
terized as a two-port network mounted in a test fixture. First, the S-parameters of the 
input (A) and the output (B) sides of the fixture must be evaluated before the S- 
parameters of the device itself can be deduced. If the test fixtures are known to be sym- 
metrical, then only one measurement is needed to determine the S-parameters of the two 
Siceseo! thesmixture. 

In the present case, however, unsymmetrical fixtures are chosen, requiring six dif- 
ferent measurements and calculations for both fixtures independently. To solve this 
problem, the input half of the fixture is first terminated at the device plane with an ar- 


bitrary load of reflection coefficient, ,, followed by other loads with different reflection 


Figure 3. 


od a _ — yal 


I 
{ 
| INPUT HALE : OUTPLT 


Q——_——<agfn—i lf (ee a 


| Ai2 | D12 
: | 


| | 
2 | 
\ wee 
Ik ON DEVICE vole 


NETWORK MEASUREMENT NEEAORE 
ANALYZER PLANES ANALYZER 
ME ASURE MENT MEASUREMENT 
ABoINs PLANE 


f 


| NETWOR¢ ! 


| 
Su = a 
ANALYZER 


: a, 


esee@eeeeeceeoeeseeeeoeeceaeaeoceeepeeaesaeeaee 8 s 8 Bg 





ee te 
Serpeeeses 
SBeeenesean 
@seeeeuaa 
epee aegcas 
SBese@eeoeas 
gSaeeenvases 
Baesupaeveeos 
ee ee 
*®esp eee ee 
@@e@aernraugs 
Beeewaane 
S@euesousen 
seeegmpese 





OI SM ® CONNECTOR 
MICROSTRIP LINE : 





mPa Gliie <2 


Sle ai = 


The Flow Diagram and Block Diagram of the Fixture. 
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coefficients. Then the output half of the fixture must also be considered with same 
method. 
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Figure 4. One Port De-Embedment 


For the one port case, as shown above, the reflection coefficient is the data needed 


since the impedance can be calculated from it: 


let Tap 
Lo iT, Zo (2.1) 


From the definition of the scattering matrix, the input reflection coefficient 1s: 


6; 
ee = ay 
(2.2) 
— 4Sy + QS17 
i ai, 
where, for load I,, 
a4 = T, 6, a) 
and 
by = S714, + S724 
(2.4) 


= S214, + S22(T 5) 


Thus, 


ay 


= ee De 
by = So) la Sal, (2.4a) 
Using the above, we get: 
= Sq 4 ae 2.5) 
HW (= S2F)) ( 


Here, S,, can be found by measuring a matched load in the fixture (T, = 90). If we 
assume that this measurement gives T, =I, then we have S,, =[,,. Using the resmin 
in Eq. (2.5) 


(Tin — Pap - Sel) = Sy2Sa 0, (2.6) 
From this, 
S1292 1 20h cee — nee) — ree vere) (2.7) 


Two more measurements with known loads will provide enough data to determine 
S,, and (S,,S,,). Suppose we measure with an open-circuited load (T,=+1). and a 
short-circuited load ([, = —1). Then, 


S12S91 + S22(Foe — Vy) = Toc — Taz (for T,=+l) (2.8) 
—$19S9) — Sx2(P se — Ty) =P se - Tu (for T,=-1) (2) 
These can be solved for 
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0c muse 


S39 = 


T9¢ — Fay) (aye — Use) 


Shige Sl 
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(2.11) 
With these two quantities known, the equation for de-embedding a general unknown 
can be written from Eq. (2.7) 


(T,, aii Myyz) 


(,=- _- 
tS 4989) + Sy9 (Pin — Tage) 


(2.12) 


But in the actual case, it is very difficult to construct the matched load and ideal 
shorted and open terminations required by this procedure. Therefore, we rearrange Eq. 
(2.12) to matrix form for three different known loads and their reflection coefficients. 
MEO) 4.(2.5)5 





Vin Si $1252} 
rT, a l a SyT, 
ie) 
—— S14 = SySy2 + S422) 
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l i ST, 
Rearranging this gives: 


Din Ss) 


a = ie + VinSa2 — (Sy1$22 — S12521) 
P (2.14) 
—s tr a VinS27 A 
where, 
A = $41522 — S12521- (2.15) 


This equation then involves the three unknowns S,,, S,,, and A. These three unknowns 
can be determined by making three independent measurements of I",, with three different 


known loads. In matrix form, the problem can be expressed as: 


ie 
td I, ] =| Si; inl 
Cl : a 
] Vin 
eee) S215 oer. (2.16) 
Ding 
73 Ping A ee 


Gpiindi S75. and A, solve: 
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LI LI 
a Ving 
Sy | = Pe Ding —l me (ze) 
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[,,, [,, and f,, are three different calculated reflection coellicients at the termunare 
of the three different loads. Also, T,,, [,, and 1,, are three tellection coeiiicienr. 
measured from the N.A. with respect to the three different loads. Now the problem re- 
quires accurate values of I’,,, ',,, f,, and T,,,, l,,, f,,,. We consider inthe next caaercn 


the computation of [,,, [,, and I’,, for reference microstrip line loads. 


it. MICROSTRIP LINE FOR THE FIXTURE 
A. MICROSTRIP LINE 


The definition of an effective relative dielectric constant, ¢,,, allows us to treat 
microstrip line like a TEM system. Relations between line and freespace parameters 
then hold (which are typical of TEM systems). With e¢,, replacing the conventional rel- 
ative dielectric constant, 

















a _ ae (3.14) 
ne = == (3.15) 
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ra A eee (3.1d) 


where /,, and v,, are the wavelength and velocity of the waves on the microstrip line and 
A, and c are the free space values at the frequency of operation. 

In Eqs. (3.1), L,, and C,, are the inductance and capacitance per unit length of the 
microstrip line. Z, is its characteristic impedance, and Z,,,, is the value Z, would have if 
the dielectric layer were replaced by air. To the lowest order of approximation, the value 
of ¢,, would be the average of the relative dielectric constant of the substrate and that 


of air, namelv: 





E (3.2) 


More exact values have been obtained, largely through solutions of the quasistatic 


model, corrected to match measured data. 
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Figure 5. Microstrip Line Dimensions. 
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If microstrip were a true TEM system, the velocity of waves on the line, the 
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impedance Z, and other line parameters would be independent of frequency. The 
quasistatic or zero frequency values of the parameters are good approximations at fre- 
quencies up to the order of a few gigaherz. When the wavelength becomes comparable 
with the cross sectional dimensions of the line, v,, and the line parameters are no longer 


frequency independent. Therefore, dispersion of effective relative dielectric constant ts 


epe(f) = &, aay | (3.34) 
es = 


where h is the substrate height, fis frequency and yu, 1s the permeability of free space. 


B. MICROSTRIP STUB-RESONANCE EVALUATION OF DIELECTRIC 
CONSTANT 
It was necessary to determine the substrate dielectric constant e, of the microstrip 


line substrates used for measurements in the present work. 


2 


Uiiessteps are: 


1. Measure the resonant frequency of a quarter wave length stub. (Substrate height 
is h= 1.5 mm, line width ratio u = 3.0 (where u = w/ h)), 


4 Different microstrip stubs, with: 

— Length = One quarter wave length 
— Substrate height = 1.5 mm 

— w/h (Line width ratio) = 3.0 


connector N connector 





Figure 6. Measurement Setup for Dielectric Constant Measurement 


2. Measure the value of S,, transmission coefficient, with respect to frequency using 
the N.A., 


3. Find the frequency (f,) at which maximum dB transmission loss due to resonance 
occurs, as shown in Fig. 7, 


S21 
S24 Magnitude plot 





Fs Frequency 


Figure 7. Determination of the Max. dB Loss Point. 
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4. Calculate e¢,, from Eq. (3.la). It is normally required to account for the fringing 
capacitance at an open microstrip termination. Computation by standard methods 
(Ref. 3: p. 138] showed that the open-end length correction was less than 1.0% for 
stub lengths greater than 50 mm. This had negligible effect on the calculated ¢,, 
and was neglected. 








dg = Ams Ere (3.4) 
a Z 
ye V €re ers 
=foAms 
therefore, 
C 
Lae a i (3.6) 
ON ere 
where 
Ams =H. 


From the equations above, 


oe aise (3.7) 


Through the several measurements as shown in Table 1, the average «,, ob- 
tained by calculation has the value of 2.833., 


Table 1. THE MEASUREMENT OF EFFECTIVE RELATIVE 
DIELECTRIC CONSTANT 


XX 







5. From stub resonance data, the corresponding ¢,, 1s calculated. Eq. (3.3) was in- 
verted by trial and error to find e, = 3.516 from the experimental e,,. 
C. CHARACTERISTIC IMPEDANCE OF MICROSTRIP LINE 
Expressions for the characteristic impedance of microstrip lines have been developed 
through numerous successive approximations, based principally on TEM models having 
the specific microstrip line cross section, with correction terms added to obtain agree- 


ment with experimental observation. 


zZ = in| (0.99u + 1.795)(0.99-4u + 5.80) (3.8) 
te p Ut Ie e5) 

Where the synthesis formula for linewidth ratio u = w/h is: 

vee / (0.607 +248). 1.321(1 +5) (3.9) 

with 

4m| exp(Z, er | (3210) 


From Eqs. 3.3, 3.6 and 3.8, one can easily get the characteristic impedance for each 
microstrip section. For example, for the given characteristic impedance of 50 Q , the 
microstrip line width ratio is u = 2.283 and the width is 3.43 mm. In the construction 
of microstrip stub loads, it 1s necessary to choose high characteristic impedance for ac- 
curate measurement, since high impedance leads to narrow linewidth and a well-defined 
point of attachment of the stub. Construction problems arise in making microstrip line 
with very small width. For example, if we want to make a microstrip line stub of 90 Q 
characteristic impedance, the required width is about 1.2 mm for given h = 1.5 mm. It 
is difficult to handle the adhesive foil material used for the microstrip itself. For this 
reason, a compromise value was chosen, using 2.0 mm width of microstrip stub for given 


substrate height h = 1.5 mm. Then, 
A=1.5mm 
w= 2.0mm 

Ww 


so = res eS Be 


Therefore, 
2. = 08.202, 


which was chosen as a working value for stub impedance. 
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Figure 8. Experimental Setup for Characterization of the Fixture 


D. STUB LOAD IMPEDANCE AND REFLECTION COEFFICIENT 


Consider the transmission-line impedance-transformation law: 


Z, cos Bi + jZ,. sin Bl 
L B JLos in B (3.11) 


Zin = Los Zo; cos Bl + Z, sin Bl 


From Eq. (3.11), when Z, goes to co (open ended), the input impedance of an open - 


ended stub becomes: 





jsin Bl ao 
Po Pou (are 
= —jZo; cot(Bl) 
using 
Zit 
are a3) 
— 
tl 
ec (3.14) 
Wee 
Therefore, 
Z,=—jZo, cot] =i si | 
2nfl./ Erp (3.15) 
605) a 


Once the Z,, is chosen, we can determine the load impedance of three different loads. 
For example, given Z,, = 68.2 ohms and microstrip line lengths of 4 cm, 3 cm, 2 cm, the 


load impedances are tabulated as shown Table 2. 


Table 2. THE LOAD IMPEDANCE OF REFERENCE STUB LOADS 








The three different reflection coefficients of the given arbitarv loads are: 


Zr — Zo 


= oS Bilis 
Zi + Lo a 


Py 


where i = 1, 2, 3. Therefore, from the given stub impedances and Z, = 50Q, the re- 


flection coefficients are tabulated in Table 3. 


Table 3. THE STUB REFLECTION COEFFICIENTS 








IV. APPLICATION OF THE DE-EMBEDDING PROCEDURE 
A. MEASUREMENT OF THE EMBEDDING FIXTURE 


The reflection coefficients required are measured by the HP 8409B vector network- 
analyzer. One can find the reflection coefficient, [,,, of the stub terminated fixtures by 
measurement of the scattering matrix element S,, using the network-analyzer. The out- 
puts of the N.A. represent the magnitude (dB loss) and phase angle (degrees) of S,,. 
Therefore, the outputs of the N.A. must be changed to complex notation using the 


simple modification: 


—Sirge Sida 
V;,=(10) 20 cos( arg S,,;)+j(10) 20 sin( arg S,,). (4.1) 
The outputs of the N.A. and T,, are obtained for each load on one fixture (A) as shown 
in Figs. 9-a to 9-c. When we consider the other fixture (B), all steps are exactly same 
as for fixture A except the fixture shape itself, shown in Figs. 10-a to 10-c. Tables for 


fixture A and fixture B data are available in Appendix A. 


B. S-PARAMETERS OF THE FIXTURE 

To consider the de-embedding calculation, we return to the original matrix equation. 
The success of this characterization technique depends heavily on the measurement of 
the load impedances used in the unterminating process for the fixtures. Once an input 
reflection coefficient has been measured with three different impedance standards, the 
S-parameters of the test fixture can be computed using the matrix introduced earlier. 
A word of caution must be offered about the convergence of this apparently simple 
problem. 

The matrix that must be inverted consists of elements very near unity, a situation 
often used to illustrate some of the numerical inversion methods. In this case, the 
Gauss-Jordan method with full pivoting is recommended. Pivoting improves the stabil- 
ity of the numerical solution by interchanging rows and columns in the Gaussian elim- 
ination. In the notation of Eq. (2.17), the scattering parameter S, is changed to 4, and 


B, for fixture A and B (respectively): 
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Figure 9. Measurement of the Load for Fixture A 
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Figure 10. 


Measurement of the Load for Fixture B 
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Ai, ie bee eee — | —| Vint 
Dy, Dy 
] Din? 
Azz | Ga) igre! Tis (4.2) 
Ving 
“A | | 13 Ding Fil cz 
where 
A= Ay, Ay — 412421. (4.3) 


The computed values for A,, and A,,, combined with a knowledge of A , can be used to 
determine the A4,, A, product. Note that the fixture is reciprocal, so A,, = A,, and for the 
other fixture, B,, = B,. Thus, the magnitudes of these terms can be calculated by taking 
the square root of the product term in the expression for A. 

The question of phase 1s left open, but is not essential. Since the A,,, 4,,, B,, and 
B,, terms always appear in product form, it does not matter which of the two possible 
phase angles one chooses, as long as the phase angles for A,, and B,, are approximately 
the same. All of these calculations are shown in Appendix E using the Fortran program 
shown. This program combines the measurement and calculated values, that is, the 
outputs of the N.A. are used in the calculation as a form of datafile. The resulting S- 


parameters of these fixtures are shown in Appendix B. 


C. DE-EMBEDDING PROCEDURE 
1, Background 

When a two-port network 1s connected between fixture A and B, the signal flow 
graph for the system may be drawn as shown in Fig. 11. Parameters AB and BA are 
mainly contributed by the finite isolation of the internal switches of the network ana- 
lvzer. The terms S,,, S;, S,,, and S,, are the scattering parameters of the device under 
test. The fixture parameters, A,,, Ay, 4,;, and Ay, can be determined by reflection 
measurements at port I, as discussed earlier. Parameters B,,, B,, and B,,B,, can be eval- 
uated by similar measurements at port 2. Two port transmission measurements with the 
device under test replaced by a through-line are used for evaluating parameters AB and 
BA, and products 4,,B,, and A,,8,,. When a through-line of length /, having Zero re- 
flections at both ends and a transmission 7 = e~#’, 1s connected between the two ports, 


the signal-flow graph of Fig. 11 can be simplified to the one shown 1n Fig. 12. 
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Figure 11. Signal Flow Graph for Measurements on a Two - Port Network 


The measured transmission 7, can be expressed as 


Ay By T 


T,= BA+ 
i l es PenB als 


(4.4) 
When the error parameters 4,, and B,, have been determined by the reflection measure- 
ments, two transmission measurements must be performed for the determination of BA 
and A,,B,,. The value of BA can be found by making a transmission measurement when 
no connection exists between the ports (that is T = 0). The parameter 4,,B,, can be 


determined from the measurement of a through - line using 


(T= BAI ~ 4,B,,e °°") 
Az, By i (4.5) 
"4 
where T = e~’, 18 being the propagation constant of the through-line. The error pa- 
rameters AB and 4,,B,, can be determined by making a similar measurement. 

Now that the scattering characteristics of the embedding networks have been 


measured, the scattering parameters V/,,, M.,, f,, and ,, can be corrected for the em- 


Z5 





Figure 12. Signal Flow Graph with a Through - Line 


bedding networks to obtain S,,, S,,, S,, and S,,. Elements M, represent the scattering 


parameters of both fixtures and device, namely: 


S14 q249)C1 — $22By1) + ($21)$)2 3B) 1412421) 


Maa ; (4.6) 
Ma, =B Ae SuPaay (4.7) 
Nyy By + SwEnBall= Suda) +(SySn4n aba) _ 
My = AB + 2B (4.9) 
where 

D = 1 — S492 — S39 Byy = Sy2891422By1 + S225 1472811 (4.10) 


24 


These equations must be inverted to extract the scattering parameters S,. An explicit 
solution for the scattering parameters of the device under test has been presented in Ref. 
[4]: p. 313. The solution is as follows: 


; Se ei) By Eee 
<< (4.11) 


E 
S2 = (4.12) 
ie 
itil < AEG ay pevala 
Sy So + ar (4.14) 
where 
N=(1 + A,.G)(1 + B,, A) — A,B, EF (4.15) 
eli AB) 
£| =—_- 4.16 
(4B) eae 
fl =—_——— ay 
(4,By) pase 
VLA,—-—A 
(412421) 
(f5. — By) 
H =-_=———"—+ 4.19 
(By > B>)) See 


In the present measurement, inclusion of the leakage paths BA and AB ts not necessary 
because of the one-port calibration procedure used. Therefore, Eqs. (4.16, 4.17) are: 


M, > 


kE=-eom 4.16a) 
A,B, ( 
Mp, 
ee (4.17a) 
A,B, 


Some variations of the above calibration procedure are available in Ref. [4]: p. 315. 


ZS 


The de-embedding procedure used in this work 1s based on Eys. (4.11 to 4.15) 
and Eqs. (4.16a, 4.17a, 4.18 and 4.19) which use the 9 key parameters: 


Sain Dios aye S43, ES E G, fi, iN 


and the foregoing derivations. The coinplete block diagram is shown in Fig. 14. If we 
know 4, and B8,, the scattering paraineters of fixtures A and B, and the parameters of 
both fixtures plus device, Jf,, then the scattering parameters of the device itself are ob- 
tuined by using Eqs. 4.11 to 4.14. In order to test the de-embedding procedure, consider 
the Sct Up Of ticmno- 
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Figure 13. Device and Fixture Assembly for the 51 Ohm Resistor 


It is necessary to choose a device that has known scattering parameters; then it will be 
possible to compare the output of the de-cembedding experiment with the theoretical 
parameters of device itself. The initial test object is a series-connected 51 Ohm resistor 


which has a small inductance so that the impedance of the 51 Ohm device 1s: 


Z=R+jwLl 
#51 + j22fL 


where we may assume the parasitic inductance associated with the resistor to be 
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Figure 14. The Block Diagram of the De-Embedding Procedure 
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The scattering matrix of an impedance Z connected in series with the trans- 


mission line is : 








Z 2 
~~ ESy) Sap] Z Z 
2+2 2+2 
where 
a AD 
< ge? 


and Z, is the characteristic impedance (50 Ohm) of the connected line. In this case, the 
TOUCHSTONE CAD program is useful in calculating the scattering parameters. 
These theoretical outputs are shown in Appendix C. Also the magnitude plots 


of the calculated scattering parameters of this device are shown in Fig. 15. 
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Fieure 15. Theoretical Scattering Parameters of 51 Ohm Resistor 


2. Device Embedment 


The S-parameters of both fixtures plus device are measured by the network an- 
alyzer at planes AA’ and BB’ as shown in Fig. 16. Typical measured magnitude data is 
shown in Fig. 17 with data tables in Appendix C. 
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Figure 16. Signal Flow graph for the Measurement of the (Fixture + Device) S- 


parameters 
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Figure 17. The Magnitude Plot of Both Microstrip Fixture and 51 Ohm Resistor 
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3. De-embedding the Device 

Now we are ready for calculation of the S-parameters of the device. The S- 
parameters of fixture A, fixture B and both fixtures plus device constitute the necessary 
data. Once the microstrip circuit board ts prepared, then the network analyzer provides 
the three different input reflection coefficients [,,, f,,; and |, with respect tomnreesai 
ferent calibration loads. Also computations are carried out to find the load reflection 
coefficients of these three different loads, [,,, ,, and I’,,. After following the same 
procedure for both fixture A and fixture B, we can get the respective S-parameters. 
From the matrix equation (Eq. 2.17), the three unknowns S,,, S,, and A are solved. For 
de-embedding the device itself, we need the additional S-parameters for the device plus 
fixture. These three sets of S-parameters are utilized to solve the de-embedding problem. 

Equations 4.11 to 4.14 are now ready for use. Through these equations, the 
de-embedded outputs are obtained, as shown in Fig, 18 (the table of de-embedded S- 
parameters is shown in Appendix C). In addition to the plot of the de-embedded mag- 
nitudes, the comparison of the scattering parameters S,,, $,, and S,, of the de-embedded 
output with the parameters of the embedded resistor and the theoretical values are also 


available in Fig. 19 to Fig. 21 (data tables in Appendix C). 
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Figure 18. Magnitude Plot of the De-embedded Parameters of 51 Ohm Resister 
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Figure 19. Comparison of the S11 of the 51 Ohm Resistor with the Theoretical 


Values 


34 


Jd8 
Loss 
20 


Comparison of Output $12 

16 (1). Fixtures + 51 Ohm Resistor 
(2). Theoretical of Resistor 
(3). De-embedded output 


12 
8 
: Se : e 
wee “7 IN ae 
—— (2) Sg ee Oe ~ a oa? — 
il a =" a PA es; pa 
2 eae ——_— - a ~ 7 Tomer ee 
ae Pl ee 
J | as Ss ae (3) 


) 
eee se 2606 29 =#6032¢«35«~ «38 «4.1 4.4 4.7 5.0 


Freq. (GHz) 


Figure 20. Comparison of the Si2 of the 51 Ohm Resistor with the Theoretical 


Values 
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Figure 21. Comparison of the S22 of the 51 Ohm Resistor with the Theoretical 


Values 
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VY. FURTHER DE - EMBEDDING METHODS 


A. MEASUREMENT OF DEVICE PARAMETERS USING A SYMMETRIC 
ELXTURE. 
1. Background 

As discussed earlier, if the device is constructed without'a coaxial or waveguide 
interface, de-embedding becomes more difficult. In some cases, it is possible to proceed 
using an auxiliary standard whose properties are calculated or measured at low fre- 
quencies, but, in other cases, no suitable auxiliary standard is available. In either case, 
if the fixture is symmetric, a de-embedding procedure that does not require the use of 
an auxiliary standard can be developed [Ref. 5]. 

The concept for symmetric de-embedding is illustrated in Fig. 22. Fig. 22-a 
shows a device inserted in a symmetric test fixture. We assume that the device has 
complementarv interfaces on the input and output sides so that when the device is re- 
moved from the fixture, the two halves of the fixture can be connected, as shown in Fig. 
22-b. The scattering parameters of each half of the test fixture can be extracted from the 
Scattering parameters of the whole test fixture connected, as illustrated in Fig. 22-b. 
Then the device is inserted into the test fixture, as illustrated in Fig. 22-a, and the scat- 
tering parameters are again measured. 

The method for implementing the concept outlined above is as follows. If we 
have a fixture whose T-parameters in a 50 ohm system are [7F]J, then if LHTF] re- 
presents the square root of [TF], 


eee \LATrl]= CTF. (Se) 


Now suppose we want to measure the T-parameters of a device in a 50 ohm system 
[TS]. If we separate the fixture into two symmetric parts and inserted the device to be 
measured, then for the combination we would measure a set of T-parameters LTFS ] 


given by: 
PPPS =TATFI\LTSILATF). (5.2) 
To find [TS], we would simply invert above equation to get 


PS =(HTF) CTFSICHRE] . (5.3) 
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Figure 22. Device Embedded Between Symmetric Fixture and Test Fixture 


This procedure de-embeds the device parameters [7S] from the measured parameters 


L7FS] and CTF]. The steps required to find [7S ] in a symmetric fixture are as follows: 


1. 


Measure the scattering parameters of both fixture plus device, [SFS], using the 
network analyzer, 


SS mS 
LSFS] = (5.4) 
SFS,,  SFS4, 
then transform to the T-parameters: 
—M, SFS), 
SFS>, Ss 
LTFS]= (533) 
SFS>, SFS>, 


where M, = SFS,,SFS, + SFS,,SFS,). 


. Measure the scattering parameters of both fixtures [SF]. In practice, both fixtures 


are measured with the network analyzer individually: 
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HSF,,  HSF,; 
CHSA = CEL] = CHS) =| | 











(5.6 
Si, tIS5, | 
Then, transforming to the T-paramieters: 
meds ie 
fe | = lo 
—HSF,, l 
HSE ELISE, 
mmeuem, — ii ST, Sf, + ASF,,7SE,, 
. Invert [HTF]: 
ATF, —-HTF,, 
M; M; 
ea. = (5.8) 
M, MM; 
mencny i — 1 1F,A1TP,— ATF, ATE. 
. Find device T-parameters [7S]: 
iim eet | Mes 1 Sis eed Opt) 
(TS]= (5.9) 
UNF oliiay TFS,, TFS; Hira fe, 
. Transform LTS] to S-parameters CSS]: 
Si ee 
[SSeS 
i S22 
where 
TS), - 
i Sits, (5.10) 
SSS e i 
Soe 15 Oa 
12 11 eS 
l 
ae SU: 
Sn = FE (5.12) 
ES; p 
= — Salo 
S2) aS wells} 
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2. De-Embedding of a Filter Unit 

As a test demonstration of the symmetrical de-embedment procedure, an 8436A 
bandpass filter (8 - 12 GHz) is used as the test device. All components of the network 
had type APC-7 test port connectors. The outputs [SFS] and [HSF] of the network 
analyzer are shown in Appendix D. To get the de-embedded outputs, the Fortran pro- 
gram shown in Appendix F is used. The magnitude plot of the filter embedded, the filter 
itself, and the de-embedded filter are shown in Figs. 23 to 25. Data Tables are shown 
in Appendix D. The comparison of the scattering parameters S,,, S,, and S,, of the filter 
itself, and embedded filter are shown in Figs. 26 to 28. Tabulated data is shown in Ap- 


pendix D. 
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Figure 25. The S-parameter Magnitude Plot of De-embedded Parameters of Filter 
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B. ISD (THROUGH, SHORT, DELAY) METHOD 
This procedure, developed by R. Speciale [Refs. 6, 7 ], is intended for svstem calibration 
and systematic error removal in automated S-parameter measurement. This procedure 
uses only three non-critical reference standards, called the “Through”, the “Short” and 
the “Delay”, and leads to a set of equations for which a fast. non-iterative, explicit sol- 
ution is known [Ref. 6]. 

1. The Basic Concept 

This procedure uses only three mechanically simple reference standards, known 
respectively as the “Through”, the “Short”, and the “Delay”. 

The “Through” is defined as the direct mutual connection of fixtures A and B, 
as shown in Fig. 29, when the uncalibrated S-parameter data for the device are acquired. 
The S-parameter data for “Through” are therefore acquired with the measurement lines 
to the device directly interconnected and the device removed from the signal path. 

The “Short” ts defined as a reference object presenting perfect shorts at the im- 
mediate interfaces where the unknown device is connected for measurement. 

Finally, the “Delay” is defined as any arbitrary length of 50 ohm line inserted 
between the interfaces where the unknown device is connected for measurement. Sum- 
marizing, for the “Through” and the “Delay”, the signal path is continuous between ports 
1 and 2 of the S-parameter unit, through the couplers used for connecting the object. 
A direct connection represents the “Through” and a finite arbitrary length line for the 
“Delay”, in place of the unknown object. On the contrary the signal path is blocked for 
the “Short” by perfect shorts appearing at the interfaces where the unknown device 1s 
nornially connected. 

iiieces-parameter matrices (S,], [5,9 and [—.S,] are acquired at each desired 
frequency during system calibration. These matrices describe the “Through”, the “Short” 
and the “Dela” respectively. Because of the device imperfections, these matrices do not 
represent theoretically perfect “Through”, “Short” and “Delay” elements, which are as- 
sumed to exist in the actual measurement. 

Now consider the error boxes A and B. Error box A is assumed to be connected 
on the side of port 1, while error box B is assumed to be on the side of port 2, as shown 


Mimipie.. 2 
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NETWORK - ANALYZER 





Port 1 
Device 





Figure 29. Block Diagram of the TSD Calibration Procedure 


From Fig. 29, the error boxes A and B are always assumed to be inserted at the 
two ends of the measurement 50 ohm lines connecting the device to ports | and 2 of the 
S-parameter unit. An algorithm has been derived wh ~h performs an explicit, non- 
iterative calculation of the S-parameter matrices [S, ] and [S,], which characterize the 
virtual error boxes A and B. 

2. The Equations and their Explicit solution 
As discussed earlier, the T-matrices can be multiplied simply with cascaded ele- 


ments. Therefore, consider Fig. 30, where the T-matrix for the connection shown 1s: 
UT ang) =CT4) LT yp] CT,]. (5.14) 
From the above, we get: 


[Tol = ieee Gals) 
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[Ta ] [To] 





Totnes 30, Device Embedded Setup 


where [7,] 1s the corrected T-matrix and L7,p,] is the uncalibrated T-matrix for the 
device. 

The T-parameters of fixtures plus device, [7,,J, are given by the network 
analvzer. But the problem is to obtain the T-parameters of each fixture, LT, ] and 
[7,1]. To do this, we consider the error boxes A and B as shown tn Fig. 31. 

From a formal point of view, both the “Through” and the “Delay” can be con- 
sidered as different, finite lengths of 50 ohm lines, cascaded between the error boxes A 
and B. With [7,,] and [T,,], the T-matrix derived from the S-parameter data obtained 
in these two calibrations, and [7,] and [7,1], the unknown T-parameters of the error 


boxes A and B, we can write: 
CTp,J a LIC ME Ga) i eae (5.16) 
LT pI = Pee IT. ). fool!) 


For the 50 ohm lines L, and L, we have: 


Corn 0 
(7,,J= , (5.18) 


penere py — 7 + /0|. 
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Figure 31. Error Boxes of A and B with Delay 


pe? @) 
LT,.J= ; F (Sa) 
gP2 


where p; =7, + 0;. 

Equations (5.16) and (5.17) together with the two equations for the “Short” 
characterize boxes A and B. Also, consider the situation for the “Short”, as shown in 
Fig. 32, we have 

by 4 Sia ee 


erm cS ey) 


bp e S53p “7 det Sp 


%p  1+Si3 ore 


2 
where S,,, and S,,,. are S,, and S,, of the “Short”, and S,,, S,, are S, of box A and Gex 
B (det means matrix determinant), where 1= 1, 2. 

Prom Eas. (5:16) anditseh) ie cece: 





Tig Hy 


(5.22) 
iy fH, 


Sips a 
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Figure 32. Short Situation of Two Error Boxes 





: H, : 
det S4 => S494 Ee (53224) 
4 
= ie . 
Soon = =-—— (5.23) 
ao Ty2R Ky 
AKG ; 
det Sp=Siip a (5.23a) 
3 
Now substitute Eqs. (5.22) and (5.22a) into (5.20) to obtain: 
H, F 
Bo eels ; 
<7 (5.24) 
1 
nS “Hy 
and also, 
SizaS214 = SiiaS224 — det Sy 
(5.25) 


H, 
= Si a Hy |suz 


Substituting Eqsm(s 23) candies a2 say mton( Sy eia: 


Ky 
hk S975 
SiB> 
Sy35 + rer 
é A, 


S13p521B = 51189223 — det Sz 


; Ky 
=| Srp + re Sip 


Hy, 
2 sinh p 
Ky 


os a 2 sinh p 


(5.27) 


(5.33) 


Here the quantities H, and K, are the elements of the matrices [/Z] and [K], defined 


Bia (0 iy) Lie 
[KJ aaa [aaa LTp,J 


and p is defined by: 


(5.34) 


~ | es 


- ~ 


Once the S-parameters of the error boxes A and B are known as given by Exe yen (oe 


through (5.27), use the general conversion formulas, 


; Ty ay 1 
[TF] = = 


T>, T45 S94 S>» | 


maimere det S =7S,,5.. — S,.95)- 


‘As discussed before, the following calculation to de-embed the device matrix, 
CTp] = Ear: CT aps] bal] 


can be carried out using Eqs. (5.22) through (5.27). 
3. Discussion 

This calibration procedure uses only three reference standards, which is the 
Imnimum necessary for characterizing the error boxes A and B. By using two delays of 
different length, for a total of 4 standards, a check of intrinsic consistency may be ob- 
tained. The standards used are simple enough to be inserted at the location of the ob- 
ject, even if this is an active device contained in sorre kind of measurement fixture. In 
this wav, all svstem errors are removed, including those introduced bv the fixture. 

The position of the reference plane is immaterial, as long as it is kept constant 
during calibration and measurement. Also, the given solution of the relevant network 
equations is an exact, explicit solution. Finally, the procedure is essentially independent 


of the specific system configuration. 


VI. CONCLUSION 


Various sources have been cited for the extraction of the scattering parameters of a 
device by removal of the error and contributions of the fixturing setup (Refs. 8,9,10,11]. 
In Chapters II and III, the de-embedding procedure using microstrip line within the 
given conditions was discussed. This project included various aspects such as the 
microstrip line fixturing setup and the construction and handling of small experimental 
boards. But in spite of the limitations of the early-model Network analyzer HP 8409B 
which was available for this work, the best possible results were obtained using repeated 


measurements on the same experimental board to obtain mean values. 


To extract the scattering parameters of the device 51 ohm resistor, the nine basic 
parameters as discussed earlier in Chapter IV (Eqs. 4-11 to 4-19) were used. These 
equations required the scattering parameters of both fixtures since they are not sym- 
metrical. That is, individual measurements and calculations were required for both fix- 
tures. [The matrix equation (Eq. 2-17) involved the three unknown scattering parameters 
S\,, S,, and A = (S,,S,, — S,,S,,) for fixture A and B. That meant each fixture required 
three equations for the matrix solution. Three different measurements and calculations 
of microstrip line were used to set up the above matrix entries, 1.e, three different meas- 
ured reflection coefficients from the network analyzer with respect to the three different 
microstrip loads, respectively I,,,, [,,. and [,,; and three different calculated reflection 
coefficients at the terminals of the three different loads I>), §,. and Wl aetere reqminces 
The computed value for S,, and S,, are combined with a knowledge of A which can be 
used to determine the S,,S,, product. Since S,,, and S,, for the fixtures always appear in 
product form, the magnitude of these terms can be calculated as the square root of the 
product term in the expression for A. 

The de-embedded scattering parameters of the device (SI ohm resistor) were found 
through experiment to be close to the theoretically calculated values, shown in magni- 
tude plots in Figs. 19, 20, and 21 and as tabulated in Table 4. As may be seen in Figs. 
19 to 21, the theoretically calculated scattering coefficients of the 51-Ohm resistor device 
_showed only a slight continuous variation over the entire frequency range of interest. 
The measured values of the embedded device showed, however, rapid irregular fluctu- 


ations with frequency, reflecting the probable occurrence of internal resonances within 
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parasitic discontinuities of the coupling circuitry. [t is the task of the de-embedding 
procedure to remove these effects from the measurement. The curves numbered (3) in 
the above Figures show the de-embedded measurements, in comparison with the the- 
oretical curves. As may be seen, the de-embedded results followed the expected values 
closely, showing fluctuations which are assumed to be due to measurement error, and 


network-analvzer calibration errors. 


Table 4. THE COMPARISON OF THE OUTPUT FOR THE 51 
OHM RESISTOR IN DB UNITS 


Theoretical bedded 





S,, and S,, terms for fixture A and B always appears in product form, considering 
their phase is not important. As long as the phase angles of S,, for fixture A and S,, for 
fixture B are approximately the same, one can choose either of the phase angles. Some 
confusion may arise, if the scattering parameters of both fixture A and B are specified 
independently of the final device (51 ohm resistor) parameters. In this case, the ap- 


proximate phase shift through one fixture cun be computed, or determined exper- 


ts 
‘sr 


imentally, and used to choose the proper phase of S,, for fixture A and S,, for fixture 
B. 


In Chapter V, the de-embedding procedure with svmmetrical fixtures was carried 
out, using as the device a commercial bandpass filter (8-12 GHz). Standard hardware 
couplers were added to represent the fixtures. Since the fixtures were symmetrical, the 
de-embedding procedure was relatively straightforward. As may be seen in Figs. 26 to 
28, the de-embedded scattering parameters of the filter were very close to the expectation 
Values, 1.e. of the scattering parameters of filter without fixtures. [t may also be noted 
that the measured S-parameters of the embedded filter (fixtures plus filter) were very 
simular to those of the filter alone, except of course for phase; indicative of well-matched 
coupling hardware. Therefore, in this task, the de-embedded results were found to be in 


good agreement with the expected values, as shown in Figs. 26 through 28. 
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APPENDIX A. THE OUTPUTS OF THE NETWORK ANALYZER AND 
THEIR REFLECTION COEFFICIENTS FOR FIXTURE A AND B 
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Figure 33. Nfeasurement of the Load 1 (Fixture A) 


Table 5. THE MEASUREMENT OF THE LOAD I(FIXTURE A) 
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Figure 34. Mfeasurement of the Load 2 (Fixture A) 


Table 6. THE MEASUREMENT OF THE LOAD 2 (FIXTURE A) 
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Figure 35. Nfeasurement of the Load 3 (Fixture A) 


Table 7. THE MEASUREMENT OF THE LOAD 3 (FIXTURE A) 
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Figure 36. Mfleasurement of the Load | (Fixture B) 


Table 8. THE NIEASUREMENT OF THE LOAD 1 (FIXTURE B) 
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Figure 37. Nfeasurement of the Load 2 (Fixture B) 


Table 9. THE NEEASUREMENT OF THE LOAD 2 (FIXTURE B) 
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Figure 38. NMfeasurement of the Load 3 (Fixture B) 


Table 10. THE MEASUREMENT OF THE LOAD 3 (FIXTURE B) 
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APPENDIX B. THE S-PARAMETERS OF THE FIXTURES 


Table Il. THES - PARAMETERS OF FIXTURE A 


0.062 + j0.136 0.646 - 30.289 0.219 - 30.013 
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Table 12. THES - PARAMETERS OF FIXTURE B 


a 


5000.0 O124 “GS 26 0.091 - 430.685 0.291 - 30.446 
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APPENDIX C. THE S-PARAMETERS OF THE EMBEDDED AND 
DE-EMIBEDDED 51 OHM RESISTOR 


Table 13. THE THEORETICAL S - PARAMETERS OF THE 31 
OHM RESISTOR IN DB UNITS. 
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Table 14. THES - PARAMETERS OF BOTH FIXTURE AND 51 
OHM RESISTOR 
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Table 15. THE DE - EMNBBEDDED OUTPUT OF 51 OHM RE- 
SISTOR DEVICE UNDER TEST. 
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APPENDIX D. THE S-PARAMETERS OF THE FIXTURE, FILTER 
EMBEDDED, FILTER WITHOUT FIXTURES AND FILTER 
DE-EMBEDDED 


Table 16. THES - PARAMETERS OF THE FIXTURE 
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Table 17. THES - PARAMETERS OF FILTER + FIXTURE 
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Table 18. THES - PARAMETERS OF FILTER WITHOUT FIX- 
TURES IN DB UNITS 
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Table 19. THE DE-EMBEDDED OUTPUT OF FILTER UNDER 
TEST 
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APPENDIX E. FORTRAN PROGRAMI FOR DE-EMBEDDING THE 31 
OHM RESISTOR 


DIMENSION X( 200) ,MDB11( 200), PSE11( 200) 
MDBI2( 200), PSE lZC200) 
MDB21( 200), PSEZNe200) 
MDB22(20@)), PSE2Z2( 200), MDB2220200) 

REAL MDB11,MDBI2 BZ) MDB2Z 2. Ib B 222 

MAG11,MAG12 ,MAG21 ,MAG22 
IMI1, Del 2 oe 
COMPLEX $11,512 (S21 S22 yaa) ee Bal 
Al1,A12,A21,A22 ,AA121,A12B12,A21B21 
Bll Bie B2 2 ee ee 
he BD Se ipa 2 2 hi es 2 
E,FF,G,H,N,AABB 
PI=4. *ATAN( 1. ) 


DO 10 T=1, 21 
F=18. 5°10. **8+1. 5°10. **8%FLOAT(I) 
X(I)=F/1000000. 


CALL FIX1(A11,A12,A21,A22,A12A21,AMAG12, PSEA12,F) 
CALL IMFS(M11 i aii 


CALL FIX2(B11,Bi2,B21,822, 812821, BMAGIZ Potuizer 


TMAG=AMAG 12*BMAG12 
TPSE=PSEA12+PSEB12 
XRE=TMAG*COS( TPSE) 
XIM=TMAG*S INC TPSE) 
A12B12=CMPLX(XRE, XIM) 
A21B21=A12B12 


E=M12/(A12B12) 

FF=M21/(A21B21) 

G=(M11-A11)/A12A21 
H=(M22-B22)/B12B21 

N=( 1. +A22"G)**(1. +B11*H)-A22"B11°E*FF 


Slis(G*( 1) +B sh aR ireh oe. 
S12=E/N 
S21=FF/N 
S22=(H*( 1. +A22*G) -A2Z2*E*FF ) /N 


QIM1=AIMAG(S11) 


QRE1=REAL(S$11) 
QIM2=AIMAG(S12) 
QRE2=REAL(S12) 
QIM3=AIMAG(S21) 
QRE3=REAL( S21) 
QIM4=AIMAG(S22) 
QRE4=REAL( $22) 


FIND A MAGNITUDE AND PHASE( ANGLE) 


101 
10 


21 


MAG11=CABS(S11) 

MAG12=CABS(S12) 

MAG21=CABS(S21) 

MAG22=CABS( S22) 

MDB11(1)=ABS( 20*ALOG10(MAG11)) 
MDB12(1I)=ABS(20*ALOG10(MAG12)) 
MDB21(1)=ABS( 20**ALOG10( MAG21) ) 
MDB22(1)=ABS( 20*ALOG10( MAG22) ) 
PSE11(01)=ATAN2(QIM1,QRE1)*180. /PI 
Poet )\—ATANZCOIM2, ORE2)*180. /PI 
Boece Ly-ATANZ( OIM3 ,QORE3)*180../PI 
PSE22(1)=ATAN2(QIM4 , QRE4)*180. /PI 


Nia cO, Ol enc t) MDBII 1) sMbB1I2(1) MDB21(1) ,MDB22( 1), 


il Poeiteu)srokl2 1), PSEZUCEyerok220n) 
PORMATC RONG TS CIA GF 7. 2) ) 

CONTINUE 

END 


Tevereverere se severerevereve seve ves Ye ve eves ve we ve ve se Veuve te Fe Fe Vere 
SVeROURPNESiMiss ODI ,D12,D21,D22,F) 


Ivevededstedetevevedevessdevesevevesevevedssesevedeserevederevedeveve 


COM whiz p21 D22 
PI=4. *ATAN( 1. ) 


INPUT DATA FILE IS DISK 9 
( BOTH FIXTURE + DEVICE ) 


READ(9,21)DB11,ANG11,DB21,ANG21,DB12,ANG12,DB22,ANG22 
FORMAT( 8(F7. 2)) 

V11=10. **( -DB11/20. ) 

V12=10. ***( -DB12/20. ) 

V21=10. #**( -DB21/20. ) 

V22=10. **( -DB22/20. ) 

RAD11=PI/180. *ANG11 

RAD12=P1/180. *ANG12 

RAD21=P1/180. *ANG21 

RAD22=PI/180. *ANG22 


QRE11=V11*COS(RAD11) 


QRE12=V12*COS(RAD12) 
QRE21=V21*COS( RAD21) 
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QRE22=V22*COS( RAD22) 
QIM11=V11*S INC RAD11) 
QIM12=V12*S INC RAD12) 
QIM21=V21*SIN( RAD21) 
QIM22=V22*S INC RAD22) 


DIT=CMPLACQRET1, Oi 
D12=CMPLX( QRE12 ,QIM12) 
D21=CMPLX(QRE21 (Ole) 
D22=CMPEX( GRE 22, O0N22) 


RETURN 
END 


Ssedederctets devote Teves Nevers sewers eee Te WeTe IE VC IS NS We eI ae Ie Ie Ie Te Ieee IC Coe IC IC sec wo ere 
SUBROUTINE FIXI(CS11 5812 (82 S223 ses tag bee eet) 


Sedevesevevovedtevedevedovevesevededevevevedevesete ve te desieve veto ve se Fes eee eee eae eT Cee ee ICICI 
POR FIXTURES 


COMPLEX A(3,3),HC3,3);,S@i3)5 Cs WG er eee 
COMPLEA STL Si2g ocean os 


CALL DATA1(A,B,F) 
CALL CPX(A) 


DO 10 I=1,3 
C(I)=A(1,1)*B(1,1)+A(2, 1)*B(1,2)+A(3,1)*B( 1,3) 
CONTINUE 


$11=C(¢1) 
$22=C(2) 
$$121=C( 3) 

95=8 11*S22-S8121 


SIM=AIMAG(SS) 
SRE=REAL(SS) 
SMAGSS=CABS(SS) 
PSESS=ATAN2( SIM, SRE) 
SPSE12=0. 5**PSESS 
SPSE21=SPSE12 
SMAG12=SMAGSS**0. 5 
SMAG21=SMAG12 


SI2=Cell*S22-Ssl Zr 0 eS 
$21=S12 

RETURN 

END 


Sevevevevedevevevesesvevesedeteseksesevedeieds 


SUBROUTINE DATA1(D,B,F) 


a @e ON FR EX ERX FB 4 Rewer ar EVN eve ey ed 4B ir ie oY 


COMPLEX D(3,3) ,BC1,3),RINI,RINZ RING Ree ee 
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CALL RLOAD(RL1,RL2,RL3,F) 
CALL RIN(RIN1,RIN2,RIN3) 


D@int )=1)/ RL 
D(1,2)=1. /RL2 
D(1,3)=1. /RL3 
De y=RINI 
D(2,2)=RIN2 
D(2,3)=RIN3 
D(3,1)=(-1. ,.0) 
Heo =C- 10) 
D(3,3)=(-1. ,. 0) 
B(1,1)=RIN1/RL1 
B(1,2)=RIN2/RL2 
B(1,3)=RIN3/RL3 
RETURN 

END 


Yerededesevedevedcscsevevsveseveskedeseseskesteseveseveseseseve 


SUBROUTINE RIN(CRIN1,RIN2,RIN3) 


Tedeverevedevevevevededestesevesevevevesevevevevovevedeveveve 
COMPLEX RIN1,RIN2,RIN3 


CALL F1DAT(DB111,ANG111,DB112,ANG112,DB113,ANG113) 


PI=4. *ATAN( 1. ) 
INPUT DATA FILE IS DISK8 
Cia URE ONLY) 
Pit 10 2 COB lly 20. ) 
Vit 2105 C-B112/20. ) 
V113=10. **( -DB113/20. ) 
RAD111=PI/180. *ANG111 
RAD112=PI/180. *ANG112 
RAD113=P1I/180. *ANG113 


QRE111=V111*COS(RAD111) 
QRE112=V112*COS(RAD112) 
QRE113=V113*COS(RAD113) 
QIM111=V111*SINCRAD111) 
QIM112=V112*SINCRAD112) 
QIM113=V113*SINCRAD113) 


RIN1I=CMPLX(QRE111,QIM111) 
RIN2=CMPLX(QRE112,QIM112) 
RIN3=CMPLX( QRE113 ,QIM113) 
RETURN 

END 


Vevevevevevededevedtetesekeleeteteseteteteininigieivichiisiciicekdedelcccdehseedetketetekccdtelesciededte 


Pree NE FIDAT(DB111,. ANGLLL DBL12, Moe > DB113,ANG113) 


Vevevetetevedetee: reve wees le Selection ese ae se ve ae ve! wwe ‘eve slecles eate pei tracer vaca pati este ve lowe se 


~J 
Ca 


20 


10 


INPUT DATAWEFILE tesbisk & 
(For FIXTURE A) 


READ(8,20)DB111,ANG111,DB112,ANG112,DB113,ANG113 
FORMAT(6(F7. 2)) 

RETURN 

END 


oe aweawaewee'a =e2woe awe Wowa wa a te" ad te Oe as ee at ee ee ot ee) ot eo a ey tS f, a Qa 
Vevededetededetede dele sete Se Leto eRe TREATED EERE PENCE EAR ERASER REE 


SUBROUTINE FIX2(511,812,521,522,S55S,SNAGi 2S zene 


yh tie iy toe air ir ir er er ert ir ir ier tire er er rl ete re tert rl teeter ieri wie eile tel el eee el wie te ele 


FOR FUTURE 


COMPLEX A(3,3),H(3,3),BC1,3),6¢3),003, oe soiee 
COMPLEX $11,512,521,822,SSS 


CALL DATA2(A,B,F) 
CAMTG PAD) 


DO 10 I=1,3 
C(I)=A(1,1)*B(1,1)+A(2,1)*B(1,2)+A(3,1)*B(1,3) 
CONTINUE 
$11=C(1) 
$22=C(2) 
S12 1=c@ey 
$SS=S11*$22-$$121 


SIM=AIMAG(SSS) 
SRE=REAL(SSS ) 
SMAGSS=CABS(SSS) 
PSESS=ATAN2(SIM, SRE) 
SPsil2=0. 5 Fokos 
SPSE21=SPSE12 
SMAG12=SMAGSS**0. 5 
SMAG21=SMAG12 


S12=(S11*522-S5121)**0.5 
921=512 


RETURN 
END 


rife sri srisz ise rises Gelaeloe tee rise ioe rise ears 

SUBROUTINE DATA2(D,B,F) 

ce ete ae ian Ge arate ar ae eis iris sige a 

COMPLEX D(3,3),B01. 3). RINT. RINZ, RINS RE eRe ce 
CALL RLOADCRL1 ,RL2 ,RL3,F) 


CALL R2IN(RIN1,RIN2,RIN3) 
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20 


Gil) =12/RL1 
DC. 2Z)=1./RL2 
Del sj=l2/RiL3 
D(2,1)=RINi 
D(2Z,2)=RIN2 
D(2,3)=RIN3 
BCs.) =(-1. sso) 
D(3,2)=(-1. ,. 0) 
Diese) =(-1.42.0) 
Ot tee) —=RINI/ REIL 
Bile 2 )=RINZ/RL2 
BC 1,3)=RIN3/RL3 
RETURN 

END 


YVeievedvereloleleleletolelellelelesete: sata ele Yeves eve ve eves eve vs eveve 


Oe aresa ACA So RINZSRENS) 


vesle we ss ves fee we s'¢ s'es' ve ewe ve 4 Je ve we ve Seong: 


COMPLEX RIN1,RIN2,RIN3 
CALL F2DAT(DB111,ANG111,DB112,ANG112,DB113,ANG113) 
PI=4. *ATAN(1. ) 


Woe l= To. --( -DBi11/ 20. ) 
Vii =O. =DB112/ 20. ) 
Vere iO. *( =DBI13/ 20.) 
RAD111=PI/180. *ANG111 
RAD112=P1I/180. *ANG112 
RAD113=PI/180. *ANG113 


QRE111=V111*COS(RADi11) 
QRE112=V112*COS( RAD1I12) 
QRE113=V113*COS( RAD113) 
QIM111=Vil1*SIN(CRAD111) 
QIM112=V112*SINCRAD112) 
QIM113=V113*SIN(RAD1i13) 


RIN1=CMPLX(QRE111,QIM111) 
RIN2=CMPLX(QRE112 ,QIM112) 
RENS=CuPLX(QRE113,Q1M113) 
RETURN 

END 


tevekedekvededetetettokietletcigigttvinitoqhdvivieicivdccdvacecdckicfccicklecldlelotcdvictict 


pe Le Leavastae: eee Ls nosis feaniee a eae 


INPUT DATA FILE IS DISK 10 
(For FIXTURE B) 


READCLO, 20)DB111,ANG111,DB112,ANGi12,DB113,ANG113 
Rename Gr 7. 2)! ) 


Uy 


RETURN 
END 


dedevedevevededededlsveveveve tevey 
SUBROUTINE CPX(A) 


Yevedevesevevevevedtetesicsevevevesic 


PEGE Reeeh ees 
COMPLEX A(3,3),W(3),DET( 2) 


CALL CGEFA(A,3,3,IP, INFO) 


CALL CGEDICA, 2.3303 bere w 01) 


Jevededevedevedededetevetedededededesededveteteiolededetevetede 


SUBROUTINE RLOAD(RL1,RL2,RL3,F) 


ale slevevlesevesteve Vales loutsalantodawanatan'suiawon es seveveveyeve alee, 
FV ED ED ED FR CD EDL ED CXL FRED ED ED rier is CDV ED EVER EL EY EVER EDL EDd EL ED FL CD FD riy 


COMPLEX RL1,RL2,RL3,Z0,2L1,2L2,2L3 
REAL Diy iZ aks 

PI=4. *ATAN( 1. ) 

C=3, *10"*""8 

ZOS=68. 2 

BRE=2Z. 635 

ERES=(ERE)**0. 5 

L1=0. 04 

L2=0. 03 

L3=0. 02 


A1=COS( 2. *PI*F*L1*ERES/C) 
B1=SIN( 2. *PI*F*L1*ERES/C) 
A2=COS( 2. *PI*F*L2*ERES/C) 
B2=SIN( 2. *PI*F*L2*ERES/C) 
A3=COS( 2. *PI*F**L3**ERES/C) 
B3=SIN( 2. *PI*F*L3*ERES/C) 
R1=-Z0S*(A1/B1) 
R2=-ZOS*(A2/B2) 
R3=-Z0S*(A3/B3) 
ZL1=CMPLX(0. ,R1) 
ZL2=CMPLX(0. ,R2) 
ZL3=CMPLX(0. ,R3) 
ZO=CMPLX(50. ,0. ) 
RL1=(ZL1-Z0)/(ZL1+Z0) 
RL2=( ZL2-Z0)/( ZL2+Z0O) 
RL3=(ZL3-Z0) /( ZL3+Z0) 
RETURN 

END 


Jeveveveveteledvedesededeteshcedevedledeiedededtetesedkestevede 
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CON FE RFANNEFWEUNUI EW WW WW WH 


THE INPUT DATA FILE FOR DISK 8 
3 DIFFERENT LOAD FOR FIXTURE A 


Tededcdenevevenevevedededededevesevevevedevedevevevevetevede se 
07 Dy) ea | ZaoS Sere 2:14 175 
62 =635.1 2.89 jig aes, ti, J eg Ge 
oe os 0 S27. =49 22 257 55 
46 126.4 fe =138 20 Ses, =5 
2 60155 3. 75. Ome 3.86 <-74 
33 lee? Seen 65.1 4.44 -156 
04 -69.9 3. 40 4.9 4.42 119 
oo = oe ego 3. 605 =54.2 Si Sis) 46 
12 Ore si6 4.55 -119.0 Sao9 =) 4 
32 ey) 5-78 15946 3.89 -66 
650 250 O27 6 ic 4,15 -120 
o7 = ieee D760 =2525 5202 175 
Sum, >. 9 5.70 -89.0 8.22 99 
73 A0425 4769 =134.7 8.295 2) 
SO 0M= 17/1 Bd Wee on 7 6.37 -79 
POe=t2ts 4.48 117.4 4.35 -147 
99 170.4 6. 84 47.9 3.49 157 
M2 6114. 2 8.48 61.4 Bao7 106 
88 Sys 2 1.57 =i a8 Sig SHO, 58 
Sit aio a o2 1158 Dis 1 
42 -137.9 Bin SIS) 62.0 6559 =82 


ric ieraus elton: Gr pict ache eter Gcoilok mir oriecia Lorto bier ee tor fore 
THE INPUT DATA FILE FOR DISK 9 
BOtn tt PSEURE AN 51 OHM RES. 


REL OR ER EN ER ER ER ER ER ER OR ER CR ER GRAN GN ON GR OR GR CR OW ON ON GO OR a 


14 =O ee 400 BO 4.40 51 
Some ol. 2 4.71 ~aliee 4.64 0 
1S 154. 3 Geo W aD So. 9 OO as 8) 
47 143.1 Dee oes LOO. 1 Doe. = 103 
i Coe? DpoemeLod. 2 Soi! Sie 
8S. 41.0 Om Zoe Uo9. 0 6.30 160 

1 a523 oy ele) geyarge| Goes Lis 
See L328 7.49 66. 6 7.54 68 
Joma O2 0 PEP AS, 19.4 Veoe 22 
92 loa. Soo - 20.9 7,04" =25 
27m 2. Deron -/9..5 7 Me S73 
O01 LZ. 6 tee 3. 3 6aao) 12) 
59m Ol. > oom 79, 3 Too L7G 
ins) oie Pay O00 135. 8 6457 +155 
27 merce? 6.56 85.4 Cos 83 
eyo) > Ihe ee, oi CAS 34.5 6.00 34 
Ome 6. 3 Sites) 2 Us: DS. Ole = 14 
» 34 Cv? 5. 50 ses a0 04 
a1 5155 eo elo 5.34 9116 
. 46 23. DOO -al On) So Oo LL 
She 44.3 Or 054 131.4 6.10 134 
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OF OF WP UUWNOWOINNYAWOOWrH FUKE 


COMM NWHAW FWO Lk, NM UWNNM & 
md 
NEOUADNHE HEHE RUUAAH 


~45 


27 
05 
86 
Si 
79 
68 
48 
43 
94 


LS 


18 
ou 


5 8. 
SOF 
2/8 
Se 
159 
. 00 
oz 
seo]. 


WNWOW Or OYINMNW A™ OO WN AEF LO 


WINE WNDNM DW FUMNW WP hn NFr Fr 


Jere ds seve ted deakedevevevese steak oka veaeote sess eve ve ve stesie ake oe 


THE INPUT DATA FILE FOR DISK 10 
3 DIFFERENT LCAD bOR Bi URE SB 


sevesededestesocesleste veuercsesesestesevovedeasss te vereokeve seven 


wa 


89 
00 
27 
40 
62 
04 
87 
os 
41 
eal 
37) 


2210 


55 


SNe) 


Jo 


Bie) S 


6) 


pow 


8 


Ou 


a a Sy ae 
=i. O- 
oO. 
euler 
34. 
a. 
0) 
S150; 
ee. 
64. 
ele 
aley fi 
=O: 
= lisish 
147. 
Thaye 
Soe 
-43. 
aa 
143. 
OV, 


OA IY HW OWUODHODAWWWOWWONND OO 


WMHNMNMNMWHEUMNFWNNNN NM HKD PLD Yo 


a2 
46 
44 
edt 


~41 
eS, 
Sos 
se) 
rou 
OW 
“eo 


ZF, 
24 


756 
.49 
49 
770 
. 43 
90 
ae) 
m2 


= 208 
=917 
= 558 
147. 


NDNNHUUWOOWONM WE NF WO UTO OF WW UI WO 10 


ok 


vial 


Bie: 


48 


sO 
09 
765 
a) 
ae 
. 68 
noo 


93 
ak 
02 


42 
on 
ao 


70 
80 
89 
72 


SO 


ROU OUANOFRFREFFWONWFUOWNEH 1 


APPENDIX F. FORTRAN PROGRAM FOR DE-EMBEDDING THE 
FILTER 


DIMENSION X(21),MDB11(21),PSE11(21) 
MDB12(21),PSE12(21) 
MDB21(21),PSE21(21) 
MDB22( 21), PSE22(21) 

REAL MDB11,MDB12,MDB21,MDB22 
MAG11,MAG12,MAG21,MAG22 
IM11,IM12,1M21,1IM22 

COMPLEX $11,$12,S$21,S22 
Ai Aden A202 
B11,B12,B21,B22,W 
Di Die, wet p22 
@imeei2 C21.C22 
Fill, Pigameneno2 
M11,M12,M21,M22 
H11,H12,H21,H22 


PI=4, *ATAN(1. ) 
DORTO  f=1,21 


INPUT DATA FILE IS DISK8 
(FIXTURE 1 ONLY) 


READ(8,20)DB11,DB12,DB21,DB22,ANG11,ANG12,ANG21,ANG22 
PORMATC@sCh 77 2))) 


V11=10. **(-DB11/20. ) 
V12=10. ***(-DB12/20. ) 
V21=10. **(-DB21/20. ) 
V22=10. ***( -DB22/20. ) 
RAD11=PI/180. *ANG11 
RAD12=PI/180. *ANG12 
RAD21=P1/180. *ANG21 
RAD22=PI/180. *ANG22 


QRE11=V11*COS( RAD11) 
QRE12=V12*COS( RAD12) 
QRE21=V21*COS(RAD21) 
QRE22=V22*COS( RAD22 ) 
QIM11=V11*SINCRAD11) 
QIM12=V12*SINCRAD12) 
QIM21=V21*S INC RAD21) 
QIM22=V22*SINCRAD22) 


A11=CMPLX( QRE11,QIM11) 
A12=CMPLX( QRE12,QIM12) 


8 i 


A21=CMPLX(QRE21,QIM21) 
A22=CMPLX(QRE22 ,QIM22) 


"S" PARAMETERS TO ''T' PARAMETERS 


B11=-(A11*A22-A12*A21)/A21 
B12= Al1/A21 

B21=-A22/A21 

B22=1/A21 


INVERSE OF FIXTURE 1 "T'' PARAMETER 


= 1/(B11*B22-B21*B12) 
C11=W*B22 

C12=-W"B21 

C21= -wWB12 

C22= W*B11 


CALL OF 'T'PARAMETERS OF BOTH FIXTURE + FILTER 
CALE IMES( Di, D12 bea 2 
(INVERSE OF FIXTURE 1) * (BOTH FIXTURE+DEVICE) 


H11=C11*D11+C21*D12 
H12=C12*D11+C22*D12 
R210 ll D2 ie D2 
B22=€ 1 2*D21te22*D22 


CALL OF INVERSE 'T'PARAMETER OF FIXTURE 2 
CALL. PIX2¢F 11,5 1256215822) 
CASCADE OF FIXTURE 1, FILTER AND FIXTURE 2 


M11=H11*F11+H21*F12 
M12=H12*F11+H22*F12 
ds Oo ad 2 9 A a 
M2 2=h 2 2 Ze 


"T'' PARAMETERS TO ''S''PARAMETERS 


S11=M12/M22 

S12= M11-(M12*M21)/M22 
S21 ia 

S22= -M21/M22 


QIM1=AIMAG(S11) 
QRE1=REAL(S11) 
QIM2=AIMAG(S12) 
QRE2=REAL( S12) 
QIM3=AIMAG( S21) 
QRE3=REAL( S21) 
QIM4=AIMAG( S22) 
QRE4=REAL( S822) 


Ck 
tw 


C DE-EMBEDDED OUTPUT 


MAG11=CABS($11) 

MAG12=CABS(S12) 

eGo i=C88s(52 1) 

MAG22=CABS($22) 
“MDB11(1)=ABS(20**ALOG10(MAG11)) 
MDB12( I )=ABS( 20*ALOG10(MAG12)) 
4DB21(I)=ABS(20**ALOG10(MAG21)) 
MDB22(I)=ABS( 20**ALOG10(MAG22)) 
PSE11(1)=ATAN2(QIM1,QRE1)*180. /PI 
PSE12(I)=ATAN2(QIM2 ,QRE2)*180. /PI 
PSE21(I)=ATAN2(QIM3,QRE3)*160. /PI 
PSE22(1I)=ATAN2(QIM4,QRE4)*180. /PI 


=5680 
XC I)=N+320. *FLOAT(I) 
WRITE(6,100) XCI),MDB11(1),MDB12(I) ,MDB21(1),MDB22(I), 
1 Poli h@@m Port 201), PSb2i(1),PSk22(1) 
MOOMNEEORMATC ° ,2X,F6.0,8(1K,F7.2)) 
10 CONTINUE 


END 

G Sedsseslcvevedevevdeseceseskeslccdeste dese ses seve deve vesdese dese vesese sess veseve ste 
SUBROUPINE. BMPs) (Dt. Di2,D21,D22 51) 

G Sevescvesesicvedevesevesesedevevedesesleseskesestesesicvecdcsdeslesesesviedeseskedese sede 


COMPLEX Al1,A12,A21,A22 
COMPLEX D11,D12,D21,D22 
PI=4. *ATAN(1. ) 


TNPURSO ATA hie los PISsk 9 
CEOTH PIxtURE + DEVICE ) 


READ(9 ,21)DB11,DB12 ,DB21,DB22 ,ANG11,ANG12,ANG21,ANG22 
ou FORMAT(8(F7. 2) ) 


Vil=10, -C-)e 11/20, } 
Vi2=l057-@-De12/ 20. ) 
V21=10. *%*( -DB21/20. ) 
Vi22=10, (=D B22/ 20. ) 
RAD11=PI/180. *ANG11 
RAD12=P1I/180. *ANG12 
RAD21=PI/180. *ANG21 
RAD22=PI/180. *ANG22 


QRE11=V11*COS(RAD11) 
QRE12=V12*COS(RAD12) 
QRE21=V21*COS(RAD21) 
QRE22=V22*COS(RAD22) 
QIM11=V1i1*SIN(CRAD11) 
QIM12=V12*SINCRAD12) 
QIM21=V21*SIN(CRAD21) 
QIM22=V22*S INCRAD22) 
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20 


A11=CMPLX(QRE11,QIM11) 
A12=CMPLX( QRE12,QIM12) 
A21=CMPLX(QRE21,QIM21) 
A22=CMPLX(QRE22,QIM22) 


D11=-(A11*A22-A12*A21)/A21 
Di2= A11i/A21 

D21=-A22/A21 

D22=1/A21 

REVURN 

END 


Jesestecevcstcledeseskevevesesedevevesedesdedesesk teakseskcsesesesese 
SUBROUTINE FIX2Z(C Eile 2 2) 


Setevcokeve de veseseveseveseseseseseseyesesesesescsese seat te aCe 


COMPLEX A11,A12,A21,A22 
COMPLEX Bil biz BZA bce 
COMBLEA Cil ,Cl23G 7 G22 
COMNPEEXN Fille her eee 
PI=4, *ATAN( 1. ) 


INPUT DATA FILE IS DISK 10 
(FIXTURE 2 ONLY) 


READ(10,20)DB11,DB12 | DB21, D822, ANGI ISANG 12, ANG2 I anGz2 


FORMAT( 8(F7. 2)) 


V11=10. **(-DB11/20. ) 
V12=10. **( -DB12/20. ) 
V21=10. **( -DB21/20. ) 
V22=10. **( -DB22/20. ) 
RAD11=PI/180. *ANG11 
RAD12=PI/180. *ANG12 
RAD21=P1I/180. *ANG21 
RAD22=PI/180. *ANG22 


QRE11=V11**COS(RAD11) 
QRE12=V12*COS(RAD12) 
QRE21=V21*COS(RAD21) 
QRE22=V22*COS(RAD22) 
QIM11=V11*SIN(RAD11) 
QIM12=V12*SIN(RAD12) 
QIM21=V21*SIN(RAD21) 
QIM22=V22*SIN(RAD22) 


A11=CMPLX(QRE11,QIM11) 
A12=CMPLX(QRE12,QIM12) 
A21=CMPLX(QRE21,QIM21) 
A22=CMPLX(QRE22 ,QIM22) 


B11=-(A11%*A22-A12*A21)/A21 
B12= Al1/A21 
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Q 


Bez="AZ2/ Ae | 
B22=1/A21 


W= 1/(B11*B22-B21*B12) 
F11=W*B22 

F12=-W*B21 

F21= -W*B12 

F22= W*B1l1 

Re cURA 

END 


Jerededesesevesevesteveskesksdslcseseskescslestslcutcle sess veseste st 
The Input Data File for Disk8 


seveveseslesescccveskevevcvedkeskevledesesesesesescckvsescvesesese 


48 sls Oro Ome 6... 72a 1 6240 
83 720 Oy Ties. 0. 54.9 
ual es O; Dimes. 40 O46. 2 
81 0) 0.14 19.94 142.8 
41 A, OP ome. oes 7 1.7 
49 09 0 Wiel, 7a 133. 1 
80 05 0 Oh ee re -) 80.8 
oi PEO Os Ome. 30179, 7 
Al ks, Opie so 5 160.7 
60 15 O.11 18.74 144.4 
23 m0 OOo. 76 shele cs: 
56 Ae. O720 16-14 54.1 
00 We Oe)  2Zorol =19 20 
90 Lo Oi 327 0 lea Ge7 
20 PAS O20 «6-26.54 °=103.0 
25 Hols, ee ae. Oo / ~ =97. 6 
41 0. Oecome). §6.-147..4 
90 | Oeeomec |. 55 =1609.6 
80 25 O: 2025.57 5 =96.5 
28 ae De 20m Opera ot. 
14 BS, Deon 0. 70e=102.3 


WODWND FR WO WIRD OS OW W ADO W OO fo 


' 
—— 
re 
i 


OPM FMM DWNW FMNONWWONNMWUMW 


COWFrRNMUOUNWAMWMWA FO OWWNHNM EF OO 


@) 


serevevesevevevedvesese ves 


feta 


sereve verde se desteve dove ve ve ve ve ve ve 


The Input Data File for Disk9 


Ievedcdvovvcdevevesevesevevedevedevededevesevevevevcvclevek 


.84 87.81 88.98 17 2 eal 0). 6 
.86 88.94 90.45 ek 2 2)5 US 
> DOmee Oreo ae 7 oe O 1.9 0ien 130.56 
af Seer Os of leo 1.03 -83.1 - 
~Jemeoe. 79 52597) Oz 45.7 
1 Seem bs O 6 Seas DOs oO. 7 
250 . 86 Peele) 2 Cash: 
2, 77 oO 81.” 16747 Sel 220 
ks) 1. 40 1.45) AAS J16.8 
roe 220 .28 (235754 81.8 - 
BLS: ee 10, Les o) loss 2.4 
203 ile e 8, 1.3040 47.5 
25 0 1.26 -12235 eco. 
65 1 OZ 97 Oe oes. 1 t= 
.43 96 100° ee Oe 23 
54 80 265 Jee06 -4,/7 
os 90 796 1/7 Umea, 9 
.34 Os LiZOne 17 S23 oe 
.74 120 1. 30 eS ce IS fe: 
. 08 2:25 1.49 6.05 -44.1 
~ 94 25 36 1.92 4,47 102.1 
vekivkkkakiededevevedsedevedevesededtek ke kevetek sek ex 
The Input Data File for Disk10 
detetedededsvedsdodevodedsiedcdede ied keicvete te ves sek teve 
7O3 =. 06 §=0) 052m 1c, OZ aS. oe 
82 #205 -=@020 50020, 28 =835.5 
sc) -.03 "0,05 a1: 97 =117956 
: 26 =. 055" -OG20Seeek. 17-7 2 
16 205. > =O70Seeel2. 2670141, 1 
ai) J) =, 05-07-0201) “15:70 Ze} 
34 -0055-0204520.535 93.6 
m8) 200 0.10 24.67 73.7 ° 
Ra, . 08 On bomeelo 47 6 41.8 - 
02 yal Or 7 Cael. 90 Lis 
. 46 ~95 S210, 21, 2256.4 
a 4.65 32 20 6.64 =139.7 
0G 2621 0.67 7.73 164.4 
A) ol On 720ml. o> 13.3 
oy , Oo 0, O6re lh? 61 sO ee 
wel Oo O33 sala 0 LE liCrec 
ae Pe k®, 0 US s7Eee08 93995 5e 
238 0 OF OS salto? ly = 12958 
eo) Os 0; OZtape 28) 2 = O28 
. 48 705 “i-0. 02) 926. 394= 11626 
. 24 is O. 10 16.524=12457 
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